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ABSTRACT Microarrays containing 1 046 human cDNAs 
of unknown sequence were printed on glass with high-speed 
robotics. These LO-cm^ DNA •'chips'' were used to quantita- 
tively monitor differential expression of the cognate human 
genes using a highly sensitive two-color hybridization assav 
Array elements that displayed diHerential expression patterns 
under given experimental conditions were characterized by 
sequencing. The identification of known and novel heat shock 
and phorbol ester-regulated genes in human T cells demon- 
strates the sensitivity of the assay. Parallel gene analysis with 
microarrays provides a rapid and efficient method for large- 
scale human gene discovery. 

Biology has eniercd the genome era (1). Complete genome 
sequences for all of the model organisms and human will 
probably be available by the year 2003 (2). Torrents of human 
expri^cd sequence lags (ESTs) provide a starting poini for 
elucidaimg the funaion of tens of thousands of cognate genes 
(3). Genome analysis will provide insights into growth, devei- 
opmcni. differentiation, homeostasis, aging, and the onset of 
diseases ( 1-3). A detailed understanding of the human genome 
Will require the implementation of sophisticated methods for 
gene expression analysis and gene discovery. 

Recently, a microarray-based method for high-throuehput 
monitoring of plant gene expression was described (4)rThis 
■'chip*'-bascd approach involved using microarravs of cDNA 
clones as gene-specific hybridization targets to quantitatively 
measure expression of the corresponding plant genes (4 5) A 
two-color fluorescence labeling and deieaion scheme 'facili- 
lated sensitive differential expression analysis of different 
plant tissues (4, 5). The efficiency of this approach for studies 
m higher plants suggested the use of this method for human 
genome analysis (4-7). Here, we report the use of cDNA 
microarrays for human gene expression monitoring, bioloeical 
investigation, and gene discovery. 

MATERIALS AND METHODS 
Human cDNA Clones. The cDNA library was made with 
mRNA from human peripheral blood Ivmphocvics trans- 
formed with the Epsiein-Barr virus. Inserts >6d0 bp were 
cloned mto the lambda vector AYES-R to generate lO'-lO" 
recombinants. Bacterial transformants were obtained bv in- 
fecting £. colt strain JMI07/AKC. Colonies were picked at 
random and propagated in a 96-weII formal, and miniivsaie 
DNA was prepared by alkaline Ivsis using REAL preps 
(Oiagen, Chaiswonh, CA). Inserts were amplified bv PCR in 
a 96-well formal using primers (PAN132 5'-CCTC- 
TATACTTTAACGTCAAGG; and PAN133, S'-TTGTGTG- 
GAATTGTGAGQGG) complemcniarv to the AYES 
polylmker and containing a six-carbon amino modification 
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(Glen Research. Sterling, VA) on the 5' end. PCR products 
were purified in a 96-well format using OlAquick columns 
(Oiagen). 

Microarray Preparation. Amino-modified PCR products 
were suspended at a concentration of 0.5 mc/ml in 3x 
standard saline citrate (SSC) and arraved from 96-well micro- 
liter plates onto silylaied microscope slides (CEL Associates 
Houston) using high-speed robotics (4-7). A total of 1056 
cDNAs. representing 1046 human clones and 10 A rabidopsis 
controls, were arrayed in l.O-cm- areas. Printed arravs were 
incubated for 4 hr in a humid chamber to allow rehydration of 
the array elements and rinsed, once in 0.2rr SDS'for 1 min 
twice in H:0 for 1 min. and once for 5 min in sodium' 
borohydride solution (1.0 g of NaBH4 dissolved in 300 ml of 
PBS and 100 ml of 100% eihanol). The arravs were submerged 
in H:0 for 2 min at PS'^C, transferred quickly into Q29c SDS 

a7i«c:"' ""^^^ ^'^^ ^"^'^^ ^^^^^^ ^^^^ 

^r!!!**!*^."'*' Probes. Tissue mRNAs were purchased 
(CLONTECH). Jurkai mRNA was isolated as described b>- 
Schena et ai (4). Probes were made as described (4) with 
several modifications. The reverse transcriptase used here was 
Superscript II RNase (GIBCO), The Cy5-dCTP was 
purchased from Amersham. Each reverse transcription reac- 
tion contamed 3.0 ^g of total human rnKSA! Arabidopsis 
control mRNAs were made by in vitro transcription of cloned 
HAT4, HAT2Z and YcsAl.23 cDNAs (4. 8. 9) usinc an RNA 
Transcription Kit (Stratagene). For quantitation, the mRNAs 
were doped into the reverse transcription reaction at ratios of 
1:100.000, 1:10,000. and 1 : 1000 (wi/wi) respectively. Following 
the reverse transcription step, samples were treated with 2.5 ^ 
of ] M sodium hydroxide for 10 min at 3rc, then neutralized 
by adding 2.5 ^1 of 1 M Tris-HCI (pH 6.8) and 2.0 ^1 of 1 M 
HCl. Probe mixtures contained cDNA produas derived from 
3 of total mRNA. suspended in 5.0 mI of hybridization 
buffer (5x SSC plus 0.29c SDS). 

Hybridization and Scanning. Probes were hybridized to 
1.0-cm- microarrays under a 14 x 14 mm glass coverslip for 
6-12 hr at 60^ in a custom-built hybridization chamber (4-7). 
Arrays were washed for 5 min at room temperature (25''C) in 
low stringency wash buffer (Ix SSC/0.2% SDS). then for 10 
mm at room temperature in high stringency wash buffer (0.1 x 
SSC/029c SDS). Arrays were scanned in 0.1 x SSC using a 
fluorescence laser scanning device (4-7). fitted with a custom 
filler set (Chroma Technology. Bratileboro, VT). Accurate 
differential expression measurements (i.e., final fluorescence 
ratios ) were obtained by tak ing the average of the ratios of two 
independent hybridizations. 

Abbreviation: EST. expressed sequence tag. 

Data deposition. The sequences rcporicd in this paper have been 
deposited in the GenBank data base (accession nos. U566.S4-US6660) 
TTo whom reprint requcsu should be addressed, c-mail: schenafi- 
cingm.stanford.edu. 
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Cell Culture. Jurkat cells were grown in a tissue culture 
incubator OTC and 5Tf CO:) in RPMI medium supplemenicd 
with lOrr fetal bovine serum. 100 of streptomycin per ml. 
and 500 units of penicillin per ml. Heat shock corresponded to 
a 4.hr incubation at 43X. Phorbol ester treated cells were 
grown for 4 hr in the presence of 50 ng of phorbol 12-mvrisiate 
13-acetaie (PMA) per ml. 

RNA Blotting. Dot blots were performed as described (4). 

DNA Sequencing. Sequences were obtained using the 
PAN132 and PAM33 primers and a 373A automated se- 
quencer, according to the instructions of the manufacturer 
{Applied Biosystems). 

Computer Graphics and Informatics. Pseudocolor represen- 
tations of fluorescent images were made with National Institutes 
of Health image software (version 1.52). Software fordiffcrcniiaJ 
expression representations was purchased from Imaging Re- 
search (St. Catherine's, ON, Canada). Sequence searches were 
made to the nonredundant nucleotide data base at the National 
Center for Bioiechnolog>' Information (NCBI) using Macintosh 
BLAST software. The EST data base was accessed via the World 
Wide Web (hup:/ wwu'.ncbi.nim.n ih.gov/). 

RESULTS 

Gene Discovery and the Heat Shock Response. Microarrays 
were used to examine the heat shock response in cultured 
human T (Jurkat) cells. Control (37*0) and heat-treated 
(43T) cells were harvested and lysed. and total mRNA from 
the two cell samples was labeled by reverse transcriptase 
incorporation of fluorescein- and Cy5-dCTP. respeaively. In 
a second set of labeling reactions, the fluorescent groups were 
"swapped" such that samples from control and heat-treated 
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samples were labeled with Cy5. and f luorcscein-dCTP respcc- 
lively. Each pair of fluorescent probes was hvbridizcd to a 
10o6-elemeni microarray. The arravs were washed at hich 
strmgenc>- and scanned with a confocal laser scanning device 
to detect emission of the two fluorescent croups. 

Hybridization signals were observed to >95^f of the human 
cDNA array elements, but not to anv of the Arabidopsis 
negative controls (Fig. 1). Fluorescence iniensiiics spanned 
more than three orders of magnitude for the 1046 arrav 
elements surveyed (Fig. 1). Comparative expression analvsisof 
heat shocked versus control cells in the two experiments 
revealed 17 array elements that displaved altered fluorescence 
ratios of ^2.0-fold (Figs, 1 and Z^V Of the 17 putative 
differentially expressed genes, U were induced bv heat shock 
treatment and 6 displayed modest repression ( Figs. 1 and ^ -I ) 
To determine the idemit>' of the hcat-reculatcd ce~nes 
cDNAs corresponding to each of the 17 arrav elements were 
sequenced on the proximal and distal end. Data base searches 
revealed perfect matches for 14 of the 17 clones, and in each 
case proxmial and distal cDNA sequences mapped to the same 
gene (Table I). Of the 1046 human genes examined on the 
microarray. the five most highly induced in heat-treated cells 
were heat shock protein 90a (hsp90a). dnaJ. hsp90i3. poK-u- 
biquiiin and t-complex polv-peptide-l (icp-l) (Table 1). Three 
of the 17 clones did not match anventrv in the public database 
though one of the clones (B7) exhibited sicnificant homology' 
to an EST from Caenorhabditis elegans (Table 1). Each of the 
novel sequences (B7-B9) exhibited -2.fold induaion (Table 1 ) 
and relatively low-level expression (Table 2). 

To confirm the microarray results. mRNA levels for each of 
the genes were measured by RNA blotting. Each of the genes 
that displayed heat shock induction, including the three novel 

4>Heat Shock 
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Table 1. 



Clone 



Microarray elements corresponding to differentially expressed genes 



Row 



Column 



Ratio 



Blast ideniitv 



Bl 

b: 

B3 

B4 

B5 

B6 

87' 

BS 

B9 

BIO 

BU 

Bi: 

B13 

BI4 

B15 

B16 

B17 

BIS 

B19 

B20 

b:i 
b:: 
b:3 



24 

1 

15 
32 
17 



14 

7 

i: 
:8 

14 
20 
30 
10 
13 
7 

:i 

3 



Accession no. 



21 
31 
8 

19 
8 

31 
4 

19 

5 
8 



16 
19 
30 
26 
IS 
30 
16 



0.5 

0.5 
0.5 
0.5 
0.5 
0.5 
2.0 
2.0 

2.4 
2.4 
2.5 
2.5 
2.6 
4.0 
5.S 
6.3 
2.0 
2.1 

2.6 
3.5 
19 



CYC oxidase III 

Act in 
CYC oxidase III 
CYC oxidase III 
CYC oxidase III 
iJ-Actin 
Nover 
Novell 
Novcl^ 
Potvubiquitin 
TCP-1 

Pnlyubtquttin 

Polyubiquitm 

HSP90/J 

DnaJ homolog 

HSPQOq 

HSP9(io 

^3:-microglobulin 
Never 

f3:-mtcro2lobulin 
PGK 
NF-kBI 
PAC-l 



J01415, J01415 
NR, X00351 * 
J01415, J01415 
J01415. J0U15 
J01415, J01415 
NR. X00351 
U56653. U56654 
U56655. L' 5 6656 
U56657. U56658 
X04K03. X(UK03 
.\52K^2. X.^2Sfi: 
Ml 7597. Ml 7597 
XO4803. Xa4803 
MI 6660. Ml 6660 
D1338K. D13388 
X0727O. X07270 
M27024. X151S3 
S.^4761. M."*0683 
U56659. U56660 
S.^4761. M.WS3 
Ml 1968. L00160 
Z47744. M55643 
LI 1329. LI 1329 



. ^,f?Z!!.^T' u*^' fluorescence ratio, sequence .den..t>. and accss.on number ol cDNAs thai man.lcsied 

a differential expression pattern with probes prepared from heat shock. (61-17) or phorbol cstcr-ircaied (BlS-23) Jurkat cells 

excenVcYcix^d^^^^ T ^"""^'^ »° ^^^-""^ ^11 .enes are nuc r 

^sl^ i^Jlv f^r *"^"«^^°"^."?'>, Accession numbers reflect the highest score for proximal and distal sequence traces, 
respectively. CYC cytochrome c TCM. T-complex poKpept.de: HSP. heat shock protein: PGK. phosphoelv'era.e kinase 
?o^v(A; ""ra"' P^^-*- phosphatase of activated cells: and NR, trace no, readable due to the presence of 

'B7 is 67'7c identical lo an EST from C clcgans (D76026). 
"^No match in the public data bases. 
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Table :. Human gene expression monitored by microarray and RNA blot ana lyses 

H^prcssion level, per ID-^ mRNAs 



Qone 


Blast identitv 


Bl 


CYC osiHic^ Iff 


B2 




33 


^ * w uKioase III 


B4 


WTL oxiaase III 


B5 


\'J\. oxidase 11] 


B6 


Act in 


B7 


Novel (weakly to D76026) 


B8 


Novel 


B9 




BIO 


Polyubiquitm 


BIl 


1 wr»i 


B12 


PoWubtquitin 


B13 


Polyubiquitin 


B14 


HSP90/3 


B15 


DnaJ homolog 


B16 


HSP90a 


B17 


H5P90o 


B18 


0^microglobulin 


B19 


Novel 


B20 


^rmicroglobulin 


B21 


Phosphoglycerate kinase 


B22 


NF-KBl 


B23 


PAC-l 





Ratio 


RNA blot 


92/46 


0,^ 


100/80 


240/120 


0.5 


270/280 


36/18 


0.^ 


ND 


76/38 


0-^ 


SD 


62/31 


Oi 


ND 


180/89 


OS 


ND 


1.3/2.6 


2.0 


0.77/1.8 


2.0/4.0 


2.0 


lJ^/3.4 


0,8/1.8 


2 2 


1-2/1.8 


0.8/1.9 


2-4 


25/89 


2.3/5.5 


2.4 


7.1/27 


0.8/2.0 


2-5 


ND 


1.7/4.3 


2J> 


ND 


75/200 


2.6 


30/120 


1.0/4.0 


4.0 


1-6/13 


0.6/3J 


5.8 


3.2/29 


0.8/5.0 


6.3 


8.6/62 


1.0/2.0 


2.0 


5.4/15 


12/2.5 


2.1 


4-5/9.5 


2.7/5.9 


2.2 


ND 


2.4/6.2 


2.6 


4.7/9,2 


1.7/6.0 


3.5 


0.65/4.7 


Oi/9.5 


19 


0.21/15 



Ratio 



0.8 

1.0 

ND 

ND 

NT) 

ND 

2.3 

2.3 

IS 

3.6 

3.8 

SO 

ND 

4.0 

8,1 

9.1 

7.2 

2.8 

2S 

ND 

2.0 

7.2 

71 



determined. 

sequences, exhibited elevated mRNA levels bv dot blot analysis 
(iabie 2). In all cases, expression ratios as deiermined bv 'the 
two procedures differed by <2.foid for the genes identified in 
the heal shock experiments (Table 2). The two assays differed 
more widely in terms of assessing absolute expression levels- 
nonetheless, absolute expression as monitored on a microarray 
typically correlated with RNA blots to within a factor of five 
(Table 2). 

^Itl?!"^ Ester Signaling. To explore a signaling pathway 
disimci from the heat shock response, microarravs were used 
10 examine the cellular effects of phorbol ester treatment 
Jurkat cells were treated with phorbol ester, harvested Ivsed' 
and used as a source of mRNA Samples of mRNA from' 
untreated or phorbol ester-stimulated cells were labeled with 
reverse transcriptase. The probes were mixed, hybridized to 
microarravs. and scanned for fluorescence emission of the two 

J?/?";^"^ ^ °^ ^« ^''^y Clements displayed 

S2.0-fold elevated signals with probes from phorbol csicr- 
treated cells relative to control samples (Fig. 2B) 

To determine the identity of the phorbof ester-induced 
genes, clones corresponding to the six array elements were 
sequenced. Data base searches revealed perfect matches for 
five of the SIX sequences (Table 1). The two most highly 
induced genes were the PAC-l tyrosine phosphatase and 
nuclear factor-kappa Bl (NF-kBJ): modest activation was 
?Ti?^ . Phosphoglycerate kinase and ^^microglobulin 
(Table 1), One remaining clone (B19) did not match anv entry 
in the public data base (Table 1). B19 displayed a il-fold 
induction and. similar to the novel heat shock'genes. a reia- 
tively low absolute expression level (Tables 1 and 2). All six of 
the phorbol ester-inducible genes displayed increased steady- 
state mRNA levels Ijy RNA blotting (Table 2) PAC-l exprcs. 
sion (Fig. 1; Table 2) defined a detection limit of -1:500 000 
for the assay. - -'w.wu 

Transcript Imaging in Human Tissues. To determine 
whether microarrays could be used to monitor expression in 
human tissues, probes were prepared from human bone mar- 



row. brain prostate, and hean by labeling each mRNA sample 
with CyS-dCTTP. In a separate reaction, a control probe was 
prepared by labeling Jurkat mRNA with fluorcscein-dCTP 
■Die four Cy5.|abeled probes were each mixed with an aliquot' 
of the f luorescein.labeled control sample, and the four mu- 
tures were hybridized to separate microarravs. The arrays were 
washed and scanned for fluorescence emission, and hybrid- 
izaiion signals for each of the tissues samples wei^ normalized 

o°f Ik' tmf f ^^""^'^ expression profile for each 
01 the 1046 clones present on the array. 

Detectable expression was observed for all 15 of the heat 

shock and phorbol ester-regulated genes in the four tissue 

types examined (Fig. 3). In general, the expression level of each 

gene m Jurkat cells correlated rather closely with expression in 

the four tissues (Table 2: Fig, 3). Genes encodinc 3-aciin and 

cytochrome c oxidase, the two most hichlv expressed of the 15 

genes m Jurkat cells (Table 2), were highly expressed in bone 

marrow, bram. prostate, and heart (Fig. 14). Expression of 

cytochrome c oxidase. hsp90a. and the novel B7 sequence was 

significantly greater in heart than in the other tissues (Fig. 3). 

DISCUSSION 

Many of the heat shock genes identified in this study encode 
.rV.l'L friction either as molecular -chaperones" 
(HSP90a. HSP90^, DnaJ, TCP-l) or as mediators o^f protein 
degradation (polyubiquitin). The identification of these sc- 
qucnces is consistent with the biochemical basis of heat shock 
induction (10-15). Proteins undergo denaiuraiion at elevated 
temperatures, and those that fail to maintain proper confor- 
mation must be selcciiveiy degraded ( 10-15). It will be inter- 
esting to determine whether the three novel heat shock- 
mducible sequences (B7-B9) mediate protein folding and 
luniover or possess some other biochemical activity. Complete 
nucleotide sequence determination, conceptual translation 
expression monitoring, and biochemical analysis should pro- 
vide a detailed functional understanding of these genes 
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Fic. 3. Transcript profiles of heat shock and phorbol esier- 
regulated genes. Gene expression levels per 100.000 mRNAs U-axes) 
arc shown for 15 genes (Table 1) in human bone marrow (red), brain 
(green), prostate (blue), and heart (yellow). Genes are grouped 
according lo expression levels M-C). 
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Phorbol esier. a poieni aaivaiorof proicin kinase C (16. 17). 
induced a set of genes distinct from those involved in the heat 
shock pathway. The most highly induced gene identified in this 
study. PAC'L encodes a nuclear tyrosine kinase that may play 
a role in regulating transcription and cell c>cle progression 
(18). NF-icBl. a second phorbol esicr-inducible gene, is an 
intensively studied member of the Rel transcription factor 
family (19-21). The Rel proteins arc activated by a large 
number of stimuli, including phorbol esters. c\ tokincs. bacte- 
rial and viral pathogens, and ultraviolet light (19-21). Modest 
activation was observed for three sequences not known to be 
inducible by phorbol esters, including phosphoglycerate ki- 
nase, 0:-microgiobuiin. and a novel human gene (B19), Ex- 
tensive expression monitoring with microarrays should assist in 
understanding how each of these genes integrate into the 
highly complex phorbol ester signaling pathway. 

It is striking that four novel human" genes were discovered 
with an array of 1000 randomly chosen clones, particularly 
because the heat shock and phorbol ester signaling pathways 
have been so intensively studied (10-21). The facile discovery 
of these sequences underscores the fact that microarrays can 
be used for gene discovery in the absence of any sequence 
information. By this approach, clones are chosen at random 
from any library of interest and only those clones that display 
interesting expression patterns are sequenced and character- 
ized. This parallel assay, coupled with a modest DNA sequenc- 
ing facility, allows high-throughput human genome expression 
analysis and gene discover)'. 

Genes that are activated or repressed b\' a given stimulus 
provide functional clues to the cellular pathway involved 
(22-24). Detailed examination of these gene expression "sig- 
natures" can provide a dynamic view of the mode of action of 
a given signaling substance (22-24). Microarrays may thus 
allow rapid mechanistic examination of hormones, drugs, 
elicitors, and other small molecules; moreover, functional 
analysis of transcription factors, kinases, growth factors, cyto- 
kines, receptors, and other gene products should be possible. 
Efforts are underway to develop mRNA amplification strate- 
gies to enable probe preparation from minute tissue samples. 
This capability might allow for high-throughput patient screen- 
ing in a clinical setting. 

The current detection limit of the assay allows monitoring of 
transcripts that represent -1:500.000 (wi/wt) of the total 
mRNA. This 10-fold increase in sensitivity compared with the 
original report (4) was achieved largely by modifying the 
coupling chemistry, which reduced background fluorescence. 
The significance of this improvement is considerable in that 
approximately half the human genes identified in this study, 
including alt four novel sequences, exhibited expression levels 
below the original detection limit of 1:50.000 (4). 

The ability to detect 2-fold changes in expression was 
achieved by the use of two-color fluorescence in the labeling 
and detection schemes, digitized data collection, and custom 
software. The importance of this capability is underscored by 
the fact thai nearly all of the genes examined here exhibited 
<6-fold changes in expression. The four novel genes, which 
showed s2.2-fold activation, were probably overlooked in 
previous screens that used conventional differential expression 
techniques. It may be possible to further improve the precision 
of the microarray assay by the use of closely related fluorescent 
analogs, such as Cy3 and Cy5. m the labeling and hybridization 
reactions. 

Microarrays offer a number of advantages over other po- 
tential high-capacity approaches to expression analysis. The 
chip-based approach enables small hybridization volumes, high 
array densities, and the use of fluorescence labeling and 
detection schemes. These features provide a set of perfor- 
mance specifications that are unattainable with filter-based 
approaches (25. 26). The use of cDNA clones provides hy- 
bridization specificity that is not readily attained with oligo- 
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nucleoiidc arrays (27-30). The parallel format of the assav 
provides a simultaneous differeniial expression readout for 
>1000 genes. This contrasts with sequcncine-based methods 
which require serial data collection for expression analysis (3l' 
3-). A commercial source of cDNA microarravs would greativ 
speed the use of a chip-based approach to expression analysis' 
The availability of large numbers of ESTs (3) provides a'rich 
resource of human cDNA clones for microarravine The 
>400.000 ESTs in the public data bases representTSc^m 
subset of all human genes (3, 33). Microarravs of thousands of 
ESTs will provide a powerful analytical tool'for future human 
gene expression studies. The -100.000 genes in the human 
genome (2, 33) emphasize the need for microarravs of greater 
density. Attempts to improve microdeposition techniques are 
underway and should allow construction of arrays containing 
a compicie set of human gene targets (http://c^gm^ianford 
edu/--schena/). Microarravs of -100.000 cDNA elements 
would allow expression monitoring of the entire human ce- 
nome in a single hybridization. This capacity, coupled with 
detailed biochemical analysis of the individual gene products 
would greatly speed the functional analysis of the human 
genome. 

Wc thank S. Elledgc (selledge^ bcm.imc.edu) for the human cDNA 
library, Oiagcn representatives for help with plasmid purification, and 
A. J. Smith and colleagues at the Protein and Nucleic Acid (PAN) 
laciiity (Stanford) for olieonuclcoiide synthesis and DNA sequencinE 
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